Enhanced diffusion by reciprocal swimming 
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Purcell's scallop theorem states that swimmers deforming their shapes in a time-reversible manner 
("reciprocal" motion) cannot swim. Using numerical simulations and theoretical calculations we 
show here that in a fluctuating environment, reciprocal swimmers undergo, on time scales larger 
than that of their rotational difl'usion, diffusive dynamics with enhanced diffusivities, possibly by 
orders of magnitude, above normal translational diffusion. Reciprocal actuation does therefore lead 
to a significant advantage over non-motile behavior for small organisms such as marine bacteria. 



In addition to its importance on our macroscopic 
world, fluid mechanics plays a crucial role in many cel- 
lular processes. One example is the hydrodynamics of 
motile cells such as bacteria, spermatozoa, algae, and 
half of the microorganisms on earth P, [2] ■ Most of them 
exploit the bending or rotation of a small number of flag- 
ella (short whip-like organelles, length scale from a few 
to tens of microns) to create fluid-based locomotion Q . 
In contrast, ciliated microorganisms swim by using the 
coordinated beating of many short flagella termed cilia 
distributed along their surface Q. 

Two physical ideas govern the fluid mechanics of cell 
locomotion on small scales. The first one is the exploita- 
tion by cells of anisotropic drag-based thrust to generate 
instantaneous propulsive forces ^2i] . The second one is the 
requirement to distribute this local propulsion along the 
surfaces of organisms in a manner that does not average 
to zero over one period of cellular actuation Q. Indeed, 
on very small scales, the inertia-less equations governing 
the surrounding fluid are linear and independent of time 
(Stokes equation), and thus any actuation on the fluid 
remaining identical under a reversal of time (so-called 
"reciprocal" actuation) cannot generate any net motion. 
This is known as Purcell's scallop theorem 0, Q . 

To overcome the constraints of the scallop theorem, mi- 
croorgansims swim using wave-like deformations of their 
appendages or bodies, be it prokaryotes, eukaryotes with 
small number of flagella, or ciliates [l], 0. For defor- 
mation of synthetic swimmers, at least two degrees of 
freedom of shape change are required [1, , or further 
physical effects need to be exploited, for example those 
leading to nonlocality (hydrodynamic interactions [^), 
relaxation (actuation of flexible filaments [l3|) or nonlin- 
earity (in particular, non- Newtonian stresses [HI). 

In contrast to large organisms able to sustain direc- 
tional swimming for long periods of times, small bacteria 
quickly lose their orientation due to rotational Brownian 
motion. If a is the typical hydrodynamic size of an or- 
ganism in a fluid of viscosity r] and temperature T, this 
thermal orientation loss occurs on a typical time scale 
T rja^ /ksT, of about one second for a 1 fim bacterium 
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in water, and tens of seconds for E. coli. On time scales 
t T, the coupling between locomotion at a typical 
speed U and orientation loss [l^, [l3| leads to diffusive be- 
havior for the cells with an effective diffusivity D ^ U'^t, 
usually much larger than that due to normal Brownian 
motion. For example, dead E. coli bacteria have diffu- 
sivities of « 0.1 /im^/s while those of swimming cells are 
at least three orders of magnitude larger [l^ . This tran- 
sition from directional motion to diffusive dynamics was 
further addressed in recent work [l3 |. 

For small organisms significantly affected by Brownian 
diffusivity, we thus have the following intriguing observa- 
tion. The scallop theorem dictates how cells should de- 
form in order to undergo non-zero time-average displace- 
ments but at long times, cells always diffuse, and thus 
always display zero time- average displacement. Would it 
then be possible that similar enhanced diffusive motion 
could be obtained within the constraints of the theorem? 

In this paper we consider the fate of swimmers under- 
going reciprocal actuation in a fluctuating environment. 
Although the scallop theorem prevents swimming on av- 
erage, we show that on time scales larger than that of 
rotational diffusion, these reciprocal non-swimmers un- 
dergo diffusive motion with enhanced diffusivities, possi- 
bly by orders of magnitude, above their normal Brownian 
diffusion in translation. This result is demonstrated com- 
putationally using Brownian dynamics simulations, and 
analytically using exact calculations for the long-time ef- 
fective diffusivity of reciprocal unidirectional swimmers. 
The different regimes obtained are also captured by phys- 
ical scalings. These new results demonstrate thus that 
reciprocal actuation, useless at zero temperature, does 
in fact lead to a significant advantage over non-motile 
behavior for small organisms such as marine bacteria. 
There is thus no rms scallop theorem. 

For a first illustration of our results, we use numerical 
computations. We performed Brownian dynamics simu- 
lations [isj of a spherical swimmer (radius a — 1 /im), 
in water at T = 300 K and during a time interval of 
100 s, with results shown in Fig. [1] The instantaneous 
velocity, U, and rotation rate, fl, of the sphere satisfy 
the dynamics: Rpu ■ (U - Ugwim) = F^, Rlu • ^ = L^, 
where Uswim is the swimming speed, Rf(7 = Swrjal and 
R-La = STTi]a^l are the viscous resistances in translation 
and orientation (1 is the identity tensor), and and 
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FIG. 1: (color online) Brownian dynamics simulation of a spherical swimmer (radius a = 1 ^im), in water at T = 300 K during 
a time interval of 100 s (5 realizations are superimposed), (a): No swimming; (b): Steady swimming at speed U = 5 /im/s; 
(c) Reciprocal swimming at speed U cosujt with U — 5 ^m/s and u = 2Dj{ {Db. is the rotational diffusivity of the swimmer, 
u) — 0.33 rad/s). Case (a) is pure Brownian motion while both (b) and (c) show enhanced diffusivities. 



are, respectively, zero-mean Brownian forces and torques, 
with correlations governed by the fluctuation-dissipation 
theorem, i.e. (F^(i)F^(t')^) = 2kBTB.FU^{t - t') and 
(L^(t)L^(i')^) = 2kBTB.Ln5{t - t'). 

Simulations were performed for three different swim- 
ming behaviors; in each case five realizations are super- 
imposed in Fig. [T] In Fig. [T^, the spheres do not swim 
(Ug^vim = 0) and thus undergo pure Brownian motion. In 
Fig. [TJ), the spheres swim steadily at speed Ugwim ~Ue 
where e is a unit vector fixed to the swimmers, and 
U — b iim/s. With these parameters, the time scale for 
thermal orientation loss is on the order of r w 3 s; we are 
thus in the regime where t ^ t, and the steady swim- 
mers show diffusive behavior with a diffusion constant 
significantly larger than the Brownian one from Fig. [T^. 

Our new result is illustrated in Fig.[TJ;, where we show 
the dynamics of swimmers undergoing reciprocal motion 
with velocity Uswim = U{t)e and U{t) = Ucosut with 
[7 = 5 /im/s and uj — 2Di^ where D/j is the rotational 
diffusivity of the swimmer (w = = 0.33 rad/s). Al- 
though the swimmers display no net motion even at short 
times (by construction the swimming speed averages to 
zero over one period of actuation), it is apparent from the 
numerical results that they diffuse much faster than pure 
Brownian motion (Fig.[T^). In what follows, we use scal- 
ing arguments and theoretical calculations to rationalize 
and quantify these results. 

How can we physically account for the increase in 
swimmer diffusion? The simplest approach involves re- 
calling the dynamics of three-dimensional (3D) random 
walks [l^ [l^. If a particle at position x undergoes a 
3D random walk where steps of size £ are followed along 
random direction during time intervals 6t, then the par- 
ticle shows no average motion, (x) — 0, but undergoes 
rms spread as (x^) ~ N£'^. Since time increases as 
t N6t, we get diffusive motion with (x^) ^ Dt with 
the diffusion constant, D, scaling as ~ i'^/St. In the 
previously-understood case of steady swimming at speed 



U, the step size is the swimming speed times the time 
step, £ — USt, and the relevant time step for change of 
direction is the time scale over which the swimming direc- 
tion is lost, i.e. dt = r, leading to the well-known scaling 

D - f/v [g. 

Reciprocal non-swimmers subject to Brownian noise 
also behave as 3D random walkers, and to estimate their 
effective diffusivity, we have to consider the appearance 
of a new time scale, namely the period uj~^ of reciprocal 
actuation over which the reversal of swimming direction 
occurs. We denote by U the amplitude of the swimming 
velocity. If the period of actuation is much larger than 
the loss-of-orientation scale, i.e. ^ t, then the step 
size is expected to be limited by the orientation loss and 
scales as £ ^ Ut, leading to diffusive motion with an ex- 
pected scaling D ^ U'^t. In this low-frequency limit, the 
effective diffusion should thus show the same scaling as 
the one for steady swimmers with the velocity amplitude 
replacing the steady swimming speed. In contrast, in the 
limit where the time for reorientation is long compared 
to the period of actuation, lj^^ <^ r, then the size of the 
3D random walk step should be limited by the swimming 
amplitude, ^ ~ J7/cj while the relevant time scale for 
change of orientation remains r, leading to an expected 
high-frequency scaling for the diffusivity as D ^ U"^ /uPt. 

We now proceed to calculate exactly the effective dif- 
fusion constant for reciprocal non-swimmers in a noisy 
environment. We consider instantaneous unidirectional 
motion with speed U{t) along a direction quantified by a 
unit vector e{t) attached to the swimming frame - this 
direction is allowed to change due to rotational diffusion. 
As the swimmer is subject to noise, its position, denoted 
x(t), follows, in the absence of inertia, the dynamics 

i^{t) = U{t)e{t)+m, (1) 

where the zero-mean noise term ^ has a magnitude set 
by the fluctuation dissipation theorem {(,{t) ■ ${t')) — 
GDksT^it — t'). Here D^^t is the Brownian diffusiv- 
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ity of the non-swimming particle {D^^t = kBT/Gnrja 
for a sphere of radius a). In the absence of swimming 
([/ = 0), the swimmer displays purely Brownian motion 
and (x ■ x) « GDh^Tt in the limit t ^ oo. When U ^ 0, 
the swimmer position, Eq. ([l} , can be integrated in time 
to give 



x(<)= f U{t')e{t')dt' + f\{t')dt'. 
Jo Jo 



(2) 



The swimming direction, e, varies in time according to 
3D rotational diffusion [13] ■ We thus expect no mean 
direction, (e) =0, and an exponential loss of swimming 
direction over time as quantified by the correlation 



(3) 



with = 2Dji and Dr is the rotational diffusion coef- 
ficient for the swimmer {Dp> — ksT /Birrja^ for a sphere). 
From Eq. ([T]) wc thus first get that (x) = and as ex- 
pected, in the long-time limit, there is no net swimming. 

To quantify the effective diffusivity, we need to com- 
pute the mean square displacements. As t ^ oo, we ex- 
pect (x ■ x) « QDt in 3D, and the effective diffusion con- 
stant, D, can thus be inferred from the limit 

Z? = i lim (x-x). (4) 

Given the integration for x, Eq. we can compute 



(x-x)(t) = U{t) 



[U{t')e{t) ■ e{t') + e{t) ■ ^{t')] At' 



Jo 

Since for any times ti and ^2 we have no correlation 
(e(ti) • ^(i2)) = 0, we obtain 

(x(t) • ±{t)) = U{t) f U{t'){e{t) ■ e{t'))dt' + 3Dk^T, (6) 



(7) 



which, using Eq. ([3]), and recalling Eq. (|4]) leads to 



lim 



The effective swimmer diffusivity, Eq. ([7|), is thus given 
by the swimming velocity correlation function modulated 
by an exponential loss (for periodic swimming. Eg. 71 
should be understood as mean value over a period) [23 1. 

With our exact calculation, we can now compute the 
effective diffusivity for some simple cases. For steady 
swimming U (t) = U, Eq. ([7]) leads to 



D = Dk^T + 



(8) 



which is the classical result [13, • In the case of har- 
monic reciprocal swimming, U{t) = U cosujt, we get 



D = D 



1 U^T 



6 1 + w2r2 



(9) 
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FIG. 2: (color online) Comparison between simulations and 
theoretical predictions for a spherical swimmer. Symbols: 
Brownian dynamics simulations for the same three cases as 
in Fig. [T] (averages of 500 realizations over a time interval 
of 200 seconds). From top to bottom: steady, reciprocal and 
no swimming. Theoretical predictions shown as straight lines. 
Top (green dash-dotted line) : prediction for effective diffusion 
for steady swimming, Eq. ([8]); middle (black dashed line): pre- 
diction for diffusion by reciprocal swimming, Eq. ((9|; bottom 
(red solid line): Brownian motion. 



More generally, for periodic swimming of the form U = 
„>o exp(mwt)|, where aq is real, we obtain 



D = D 
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(10) 



which clearly displays both scalings for ujt 3> 1 and 
LUT ^ 1 discussed above. We also get from Eq. pH)) 
that we have always have D > Dk^T- For exam- 
ple, for a periodic square swimming with U(t) = —U 
during t G (— 7r/ti;,0) and then instantaneous reversal 
U{t) = +U for t e (0,7r/cj), we have Uo = U, a2p = 
and a2p+i = — 4i/7r(2p -I- 1), leading to 



D = Dk^T + 



1 



2tlo 



tanh 



2™/ 



(11) 



In Fig. [2] we show a comparison between our analyti- 
cal predictions and our computational results. We plot 
the mean square displacement of 500 realizations of the 
swimmers with the same three cases as in Fig. [T] over 
a time interval of 200s. For the three cases considered 
(no swimming, steady swimming, and reciprocal swim- 
ming), we also plot as straight lines the theoretical pre- 
diction (where D is given, respectively, by D^^t, Eq. [S] 
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and Eq. |9]). We obtain excellent quantitative agreement, 
confirming the validity of our theoretical approach. 

Biologically, our results are relevant to the dynam- 
ics of marine bacteria. Non-marine bacteria such as 
E. coli swim using a "run-and-tumble" strategy where 
straight swimming paths are followed by random re- 
orientation events [1^. As a difference, marine bacteria 
display "run-and-reverse" (or "back-and-forth" ) locomo- 
tion where high speed swimming along straight paths is 
followed by almost complete reversal of their swimming 
direction [2^ . With no bias in the characteristics of 
the paths, this is the example of a biological reciprocal 
swimmer. 

To estimate the order of magnitude of our result, let 
us consider an elongated bacterium characterized by two 
length scales, b and a ^ b. Scaling-wise, we have 
DksT ^ kBT/Tjalogia/b), Dr - ksT /rja^ \og{a/b), and 
thus the reorientation time scales as r ~ a^/Dk^T- The 
maximum enhanced diffusivity is obtained in the low fre- 
quency limit, UJT < 1. In that case, the increase of cell 
diffusivity, in a quiescent fluid environment, is given by 
D/DksT ^ U^r/DksT Pe^ where the Peclet num- 
ber is given by Pe = all / Dk^T- For blunt swimmers 
where a ^ b, even though the log terms in the diffusion 
constants disappear, the result is unchanged. For order 
one or above Peclet numbers, the diffusive behavior of 
cells is thus expected to be dominated by all swimming- 
induced terms, including the reciprocal ones. For a ten 
micron bacterium in water at room temprature, this cor- 
responds to a critical amplitude of reciprocal swimming 
of (7 w lOnm/s, less than 0.1% of the steady swim- 
ming speed of most marine bacteria [20| . For exam- 
ple, the micron-size marine bacterium Shewanella putre- 
faciens (CN32) has an average swimming speed of 100 



/im/s and run duration of about 1 second [l9|, leading 
to an expected reciprocal diffusivity of 10 /im^/s, over 
two orders of magnitude above that given by Brownian 
motion. 

Many marine bacteria are found in high-Reynolds 
number turbulent fluid environments j2T|. The frame- 
work used here remains valid provided T is interpreted 
as an effective temperature, with an equation equivalent 
to Eq. ([3]) capturing the rotational dynamics of bacteria 
in turbulent flows. The results reported in this paper 
could thus be used to describe the effective diffusion of 
marine bacteria in intermittent or turbulent flows. Our 
work could also be adapted to describe biased effective 
diffusion and chemotaxis in presence of external fields, for 
example if we allow the reciprocal swimming amplitude, 
or its frequency, to be coupled to an external chemical 
concentration. More generally, any noisy process lead- 
ing to an exponential loss of cell orientation will lead to 
enhanced diffusion for reciprocal actuation, for example 
cell-cell collisions at high density [22]. 

In summary, we have shown in this paper that recip- 
rocal swimmers, previously believed to display a useless 
form of locomotion, undergo in fact enhanced diffusion, 
possibly by orders of magnitude, over inert bodies of the 
same size. Purcell's scallop theorem, valid in the ab- 
sence of noise, can therefore not be extended in a fluc- 
tuating environment, and reciprocal (or more generally, 
unsteady) actuation can lead to significant advantages 
over non-motile behavior for small organisms. 

We thank Tim Pedley for helpful comments. Fund- 
ing by the National Science Foundation (grant CBET- 
0746285) is gratefully acknowledged. 
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